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Abstract 

We compare systematically several classes of stochastic volatility models of stock mar- 
ket fluctuations. We show that the long-time return distribution is either Gaussian 
or develops a power-law tail, while the short-time return distribution has generically a 
stretched-exponential form, but can assume also an algebraic decay, in the family of mod- 
els which we call "G ARCH" -type. The intermediate regime is found in the exponential 
Ornstein-Uhlenbeck process. We calculate also the decay of the autocorrelation function 
of volatility. 
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1. Introduction 

The modelling based on stochastic volatility 0, Q is a natural way how to implement 
non-Markovian property of certain stochastic processes into the machinery of stochastic 
differential equations, which are Markovian by definition. The idea is old, well-known 
and straightforward and amounts postulating an auxiliary process (or processes), which 
together with our process of interest form a multicomponent process, which is Markovian, 
while each of the components itself is not. The auxiliary processes may have a direct 
physical interpretation or not, depending on the situation. 

When modelling the stock-market fluctuations, the process of interest is the price 
(which is directly measurable) and the auxiliary process is the instantaneous volatility, 
or some function of it, which is, however, purely hypothetical and in principle not mea- 
surable. We can estimate its properties from the volatility measured within certain time 
interval, but nevertheless, the instantaneous volatility remains an elusive ghost. Suppos- 
ing it is subject to a process described by a stochastic differential equation, we can solve 
the equation pair for price and volatility and deduce for example the return distribution, 
autocorrelations etc. 

Since this methodology is quite general, it is of interest not only for economists but 
can be useful in various purely physical situations. One of such examples might be the 
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stochastic resonance under the influence of coloured noise. In this paper we present some 
generalisations of previously known results, stressing the generic and the special features 
of several classes of stochastic volatility models. 



2. Systematics 



Quite general ensemble of stochastic volatility models is described by a pair of stochas- 
tic differential equations (in Ito convention) with algebraic coefficients 



dX t 
dS t 



_ q7 



S?dW lt 



a (a- S t ) S? dt + g S? dW 2 



(1) 



We denote Wit and W 2t two independent Wiener processes with unit diffusion constant, 
[dWu] 2 = [dW^t] 2 = dt. In this article we disregard the possible dependence of the 
processes W\ t and W 2 t, which was investigated e. g. in @, 0, 0, S in the context of 
the leverage effect. We also omit the possibility that either W\t or W 2 t is a Levy noise 
Q, because that would lead us to a completely different area. 

The process X t describes the quantity of interest, which is observable and if we remain 
in the field of econophysics, it is typically the logarithm of the price of a commodity or 
a security. We shall call it simply the price and allow both positive and negative values 
for it. St is the auxiliary process, representing a hypothetical quantity which will be 
called volatility. Contrary to the true volatility, which is (whatever its definition might 
be) genuinely two-time quantity, like the return, the "volatility" St is an instantaneous, 
one-time stochastic variable. This does not mean that it is completely inaccessible to 
observations, but implies that its relation to observations is not a direct one and poses 
a separate problem, not to be investigated here. We consider necessary to recall these 
trivial observations, as they are often neglected in the econophysics literature. 

Some of the combinations of the exponents bear their usual names. We make an 
overview in the following table. 
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name or acronym of the model 
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Stein- Stein 
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"Ornstein-Uhlenbeck" or "OU" 





1/2 


1/2 


Heston 
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1/2 


" GARCH" 


1 
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1/2 


"geometric OU" 





3/2 


1/2 


3/2-model 



(The names in quotation marks " " are somewhat abuses of notation, because properly 
speaking they are already in use in more or less different sense.) Besides these models 
which contain only coefficients depending on S as a power, we shall investigate also the 
exp onential Ornstein-Uhlenbeck process (acronym expOU) 0, S 0, i, M, ELI E ES Q 
Il5| ] , which in its simplest form is described by the pair of equations 



dX t = e a+St dW u 

dS t = -aS t dt + g dW 2t 



(2) 
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Again, X t is the logarithm of price and St the hypothetical hidden variable. The instan- 
taneous volatility in this case is e 2 ^ rT+St \ 

If we were interested only in the auxiliary process St, we would need to solve only 
the second of the equations ([I} . Processes of this type were solved in many contexts and 
are considerably easier compared to solving the coupled pair (X t ,St). The older works 
relevant to our study are e. g. and [ItJ ■ In the former, the special case a = 0, j3 = 1, 
(7 = is solved exactly using the Lie-group technique. The latter treats approximately 
the difficult case a = 1, /3 = 1/2 which we avoid completely, as we never found it used 
in the stochastic volatility context. 

Let us make a very short historical remark at the end of this section. The study of 
stochastic volatility models in economy can be traced back to the pioneering works of 
Vasicek [l8| and Cox, Ingersoll and Ross [l9| on interest rates. In fact, the same problem 
was investigated, on the level of Fokker-Planck equation, much earlier by Feller [20j |. 
Hull and White [2l| formulated the problem on a more general basis. These results were 



adapted for the modelling of price fluctuations by Scott [9| 
and Stein 23]. The Stein-Stein model was thoroughly studied 



Hull and White 

3, 



, Stein 
although it 



exhibits some unnatural features, namely negative volatilities. These are absent in the 
Heston model [26| . which became very popular 24| 27, 28, 29|, 3(1 31, 32|. Of course, 
the potential of stochastic volatility models was used for better prediction of investment 



risks, see e. g. [H, H, 0, [H, [lj|. For further variants of the stochastic volatility 
models used in econophysics see e. g. references [34, 

M M M 0, M Hi El S3, El . 



3. Long-time properties 

The properties of the process (P) can be studied directly by the method of moments. 
This way we can easily find indirect evidence on the asymptotic distribution of price 
change. Indeed, we shall find that in the most important cases the whole set of asymptotic 
moments can be calculated and that they fall into two groups. Either the moments imply 
the Gaussian distribution of returns, or some of them diverge, indicating power-law tail 
in the return distribution. 

On the other hand, the method of moments has also serious disadvantages. Most 
importantly, it does not allow an easy implementation of the requirement of positive 
volatility. For 7=1 this problem does not occur, because the "physical" volatility is 
the square Sf , but for 7 = 1/2 we must consider it. When solving the set of Fokker- 
Planck equations, it is imposed as a boundary condition. On the level of moments 
we do not have such possibility. Therefore, we either use the fact that the process 
itself prohibits trajectories connecting positive and negative volatilities (this is the case 
of > 0) or assume that the events with negative volatilities are rare enough to be 
neglected. Typically it means that a 3> g/\/a. This is how we shall treat the case (3 = 0. 

The second disadvantage resides in the fact that we are unable to calculate the mo- 
ments for all values of the exponents a, /3, and 7. If all of them are rational, we can 
bring them to a common denominator, say, q, and calculate moments which are integer 
multiples of 1/q. For large q it is possible but difficult. However, for irrational and 
incommensurate values of the exponents this way is unfeasible. 

The general moment is defined as n m ,n{t) — (X™ &t) an d its time evolution according 



3 



to the process (HJ) is contained in the equation 



(t) = na (a fj, m<n+a -i(t) - jU m , n +a(i))- 



1 g 2 

+ -m(m - 1) /lt m _2,„ + 2 T (t) + _ 1) Mm.n-2+2,8 (*) • 



(3) 



We do not have a solution in a closed form in the general case, but we can investigate 
each model separately. We shall follow essentially the order in which the models are 
listed in the table above. 

We start with the Stein-Stein model and proceed from small values of m and n to 
higher ones. For m = the equation is 

<? 2 / 

(iO,n + na /j 0: n = naa Mo,n-i + y n ( n ~ 1) Mo : n-2 (4) 

and we can see that knowing /Uo,n-i an d /Uo,n-2 we obtain /zo,n by trivial integration. 
The starting point /zo,o(i) = (1) = 1 is obvious. The next two moments are 



Mo,i(*) = cr + (mo,i(0) - cr) e 



- a /- 



Mo,a(*) = ^ + f- + A*o,a(0) - ( a 2 + 2- 
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2 , .9 



2a 



e" 2at + (5) 



+ 2a (/i ,i(0) - <r) (e- at -e~ 2at ) 



and we can continue further as far as we please. The generic feature of the time evolution 
of moments, seen in ([3]) is an exponential relaxation to a stationary value. Setting 
Ao,n = in (01) we obtain a chain of equations for the stationary values of the moments 
|Lto, n (oo). First two of them, /xq,i(oo) = cr, ^0.2(00) = a 2 + |- can be seen in the formulae 
(O, and the general result is 



n\ (2fc)! /9 2 \* „- 2l 



confirming that the stationary distribution of St only is Gaussian with average cr and 
variance g 2 /(2a), as expected. (The fact that the distribution is Gaussian follows immedi- 
ately from the comparison with the formula for 2fc-th moment of a Gaussian distribution 
with variance s, which is [{2k)\/(2 k hi)] s k .) 

The moments with m > are just a bit more complicated. We need to investigate 
only even m, because all moments for odd m are zero. Integrating the equation ([3]) we 
can get the moments for higher and higher 777 systematically. The most interesting part 
of the result is the behaviour of the moments for large times. For example, we find, using 
©, that 

M2,o(*) = M2,o(0) + J Mo,b(t) dr ~ (a 2 + |-) t , as t -> 00 . (7) 

Generally, it is possible to show that (i) A*2/,n(i) h as a finite limit when rj — ► 00 and 
after some easy algebra we arrive at a formula summarising the long-time behaviour of 

4 



all moments 

(2l)\ / Q 2 \l 

Km r l m ,n(t) = V + 2a) ' (8) 

This result indicates that for long times the logarithmic returns have Gaussian distribu- 
tion with variance (<r 2 + |-j t. 

Analogous steps ought to be taken when we want to find the moments in the case 
of the "OU" model. Note that the moments including only the process St, i- e. those 
with m = 0, are identical to the Stein-Stein model. Therefore, the formulae ([5]) and © 
remain in force. The long-time asymptotics of the remaining moments is obtained in the 
same way as Eq. ((SJ). This time, the result is 



ljm t 1 H2i,n(t) = crVo,n(oo) (9) 



and again suggests Gaussian distribution of the returns. 

The next in line is the Heston model. The moments evolve according to the set of 
equations 

g 2 1 

jU TO ,„ + nafi m ^ n = na{<j + — (n - 1)) fi m , n -i + -m{m - 1) /Lt m _ 2 ,„+i (10) 
and solving them for the lowest non-trivial moments we get 
Mo,i(*) =(T+ (mo,i(0)-(t) e~ at 

A*a,o(t) =M2,o(0) + <rt+i(Mo,i(0)-T)(l-e- Q *) . " ' ' 

The long-time behaviour is again obtained assuming that the derivatives (^)'/i2/,n(0 
have finite limit for t — > oo and checking a posteriori that this assumption does not 
contain any contradiction. We can do even better than we could for Stein-Stein and 
"OU" models. The limits of the moments are expressed by the closed formula 



f2ZV , /o 2 \» r(^f +n) 
lim t - l fi 2ln (t) = ^-a l (^- ) V1 Z r ' (12) 



showing again without any doubt that the long-time distribution of returns is Gaussian. 

The same procedure is effective also for the "GARCH" model, where the moments 
evolve according to 

Am,™ +na(l- |^(n - l]jfi m ,n = nacr n m , n -\ + ^m(m - 1) fi m -2,n+i ■ (13) 

The evolution of the lowest moments Ho t i(t) and /i2,o(^) follows exactly the same law 
(fTTj) as in the case of the Heston model. Moreover, we have 



2a a 2 / . , 2a a 2 \ /9 „ 2, t 

M0,2 * = ~ 2+ ^0,2 - 2 ) e " ( 9 ] + 

Za — g z \ 2a — g A / 

a — g 2 



(14) 



We must be more careful here when we compute the asymptotics of the moments. Some 
of them diverge exponentially, rather than algebraically, with t, becoming effectively 
infinite. We find that the limit 



: «'.»<')=w'' + "[n('-4)]" < 15 



is finite as long as I + n < 1 + Otherwise the moment should be considered diverging. 
The interpretation of this result is, that the distribution develops power-law tails. It 
should hold both for the marginal distribution of returns, described by the moments 
with n = 0, and the marginal distribution of volatility, corresponding to the set of 
moments with I = 0. Hence, the tail exponent of the return distribution is 

r = 3 + ^. (16) 

g 2 

Therefore, the exponent is always bounded from below by r > 3, a result which is well 



consistent with the emp irical findings [44j, |45j. Note also that the empirical volatility 
distribution 46, 47, 48[ 4^, Q shows log-normal central part, followed by a power-law 
tail, which is consistent with the power-law tail in the distribution of the hypothetical 
volatility St , as is evident from (|15l) . All of that suggests that the "GARCH" member of 
the stochastic- volatility family of stock market models is close to reality. 

The next two models in the list are the "geometric OU" and the 3/2 models. Unfor- 
tunately, the method of moments does not lead to a recursive chain which would allow 
computing all the moments. The only way forward is to solve directly the correspond- 
ing Fokker-Planck equations. We shall see later how we can bypass it in the short-time 
approximation. 

The treatment of the expOU model using the method of moments is slightly more 
involved, as it requires expansion of the exponential in powers and later resummation of 
a series back, to get a compact result. The form of the equations ([2]) ensures that the 
odd moments are zero, n m ,2k-i = A*2i-i,n = 0, if the initial condition is Xq = Sq =0. In 
that case, we find, for the even moments 

,^ '"WW = f| »P[^4)] »(£)'■ (IT) 

This type of formula is already familiar to us. The interpretation is that at long times 
the logarithm of price is Gaussian-distributed with variance iexp(2<r + g 2 /a). 

4. Autocorrelations 

The autocorrelation function of returns 

C tetq (At;5t) = (\(X t+St - X t )(X t+At+St - A At )H (18) 

is a four-time quantity, not easy to calculate by the method of moments. Instead, we shall 
investigate a simpler quantity, which is supposed to share relevant properties with the 
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autocorrelation function (|18[) . Indeed, for short time intervals St the absolute increments 
of price are proportional to and we can consider the autocorrelation function 

C q (At) = (\S t S t+A t\ qy ) (19) 

as an approximation for C ro t q of St — > 0. In stationary state it is a function of one time 
argument only and we obtain for it exactly the same differential equations as we have 
had for the moments. To this end we define a slightly more general autocorrelation 

CuA*t) = (s? % +A j (20) 

and obtain the following equation for it 

d °dA^ = - va ( c u,v+ a {At) -aC u , v - 1+a (At)) + £v(v - 1) C u ^ 2+20 (At) . (21) 

We can see that the autocorrelation function with v = 1 does not depend on the value of 
f5. Especially, all the models examined in detail above have a = and the solution for 
all of them is 

C uA (At) = o-/x ,„(oo) + (w),u+i(oo) - cr/i 0)U (oo)) e" aA * . (22) 

For v > 1 the value of (3 is relevant. We could proceed iteratively finding the cor- 
relation function for v — 2, then for v = 3 etc. However, it is possible to express the 
solution as a linear combination of exponentials and then solve the set of equations for 
the coefficients. Thus we find the full solution in relatively compact form. In the case 
a = f) = (i. e. Stein-Stein and "OU" models), we get 

(2n)! ( wi 



^ 2 n n\ \ 2n 

n=0 



\2a) 



fv- m e^ v - m ^ aAt (23) 



m=0 

where /; are constants determined from the initial conditions C u-V (0) = /xo jU +i,(oo). 
Similarly, for the Heston model ((3 = 1/2) we have 

The case /? = 1 ("GARCH") is more complicated. We could also try to find general 
expressions like ((23)) and ([24|. but the decay rates of the exponentials are not ordered 



in arithmetic sequence {^ct}^ , but form a two-parametric set {ha — hg 2 /2 | lx,l% = 
0, 1, 2 . . .}. Some of these rates are negative, which indicates that some of the correlation 
functions do not exist in stationary state. This feature stems from the divergence of 
higher moments of the volatility in stationary regime. However, the general conclusion 
holds, that if the stationary autocorrelation function does exist, it decays as a combina- 
tion of exponentials. 

We can conclude that the common feature of all stochastic volatility models inves- 
tigated here is the exponential decay of autocorrelations. This is natural, because the 
process St has Markov property. For practical use, we note that the multiple timescales 
found in the expOU model Uj, [llj can successfully reproduce the slow decay of auto- 



correlations found in empirical data. Special cases when the autocorrelation does decay 
as a power law are examined e. g. in [50], but these cases do not seem to have much 
relevance for stock-market modelling. 



5. Short-time properties 



From the point of view of return distributions, the long-time properties of the stochas- 
tic volatility models discussed so far do not exhibit much diversity. The return distribu- 
tion is mostly Gaussian, with the exception of the "GARCH" model, where power-law 
tails develop. At short times, however, more variable results are found. Indeed, the em- 
pirical return distribution tend all to a Gaussian for large enough times, so it is the range 
of medium time distances at which the interesting fat tails are observed. So, it is well 
sensible to discuss the special features of the stochastic volatility models appearing at 
short times, i. e. at times shorter than those at which Gaussian limit is reached. We have 
seen that the typical decay time of autocorrelations is 1/a. We expect that this is also the 
time when the Gaussian distribution starts to prevail. Therefore, we shall investigate the 
return distributions for time distances At a -1 . This kind of approach is sometimes 
related (5l| to the Born-Oppenheimer approximation used in quantum chemistry. 

5.1. Volatility distribution 

The ingredient needed for these investigations will be the stationary distribution 
for the volatility process St- It is easy to find it, because the instantaneous volatility 
influences the price process X t , but is not affected by it, so we can find the distribution 
for St separately. We must solve the Fokker-Planck equation 



^-Ps,t(s)=a-l-[(s-a)s a P s , t (s)] + 
dt os 



in the limit t — ¥ oo. Contrary to the method of moments, we can obtain the solution 
generally for any value of the parameters a and /3. It is also easy to implement the 
reflective boundary condition at St — 0, enforcing the positivity of the volatility simply 
by multiplying the solution of f|25[) by the Heaviside function 9(s) and recomputing the 
normalisation factor so that J °° Ps(s)ds — 1. This will work for all solutions with zero 
total current. 

We can get easily the generic solution with zero current in the form 

P s (s )=M6(s ) 8 -W exp (- — — + 26 

V g z 2 + a — 2p g z I + a — 2/3/ 

where Af is the appropriate normalisation constant. The Gaussian distribution is included 
here, for the combination of parameters giving a — 2/3 = 0. 

There are two exceptional cases, marked by the vanishing denominators in the expo- 
nent in Eq. (|21)|) . When 1 + a — 2(3 = we have 

Ps(s) ^NO^s- 2 ^ 2 ^'^ exp(-^s) . (27) 

We shall call this case "Heston type", as the Heston model belongs to it. The second 
exceptional case occurs when 2 + a — 2/3 = 0. The probability density for the volatility 
is then 

P S (S) = M6(S) S -^-2a/9 2 exp ( - ^) . (28) 



This case will be named "GARCH-typc" , because it includes the "GARCH" model stud- 
ied above. 



5.2. Return distribution 

We can neglect fluctuations of the volatility during time delays At <C or 1 and con- 
sider the volatility fixed with distribution given by (|2"o| . or, in special cases, by (pM)) or 
(f2"T)) . Thus, our starting formula will be 



P A x(A,;At) = ^exp(^)p s (. 



ds 



(29) 



Contrary to the method of moments, we can proceed directly with the calculation for 
arbitrary values of the parameters a, /3, 7. We shall explore this freedom for a and /3, 
bit for the parameter 7 we use only the most studied values 1/2 and 1. The general 
strategy will be to use a suitable change of variables to perform the integral (|29p by the 
saddle-point method. 



5.2.1. 7 = 1/2 

In the general case (|2l)]) we define 



i 



(Ax) 2 



2+0 



3+0 



77 = a 



G 



At 

■ 2 (Ax) 2 



4c 



£1 

4a 
2 



5 2 (2 + < 

3+n-20 (2 



2/3) 



2/3) ? 



1-4/3 

and express the integral (|2"9"|) as 



At 

V 

4a At 



1 + a - 2/3 



(Ax) 



2/3) 



1-40 
2(3 + c,-20) 



Pax (Ax; At) 



N G 
V2nAt 



exp 



du 



where 



2 2(2 + c-2ff) 

'1-4j8 -)- u T=33 



1-40 



(30) 



(31) 



(32) 



il)(u;rf) = u I - 4 f + u 1_4p —rju 

The saddle point approximation is appropriate for large returns, |Ax|, which is just the 
regime of interest. The location of the saddle point cannot be found exactly, but it can 
be expressed as a series in powers of 77. But as we can see, 77 decreases to zero when 
Ax — > 00, so the tails of the return distribution are found both in the limit £ — > 00 and 
i] — > 0. Systematic loop expansion then gives the correction of the order 0(l/£). Finally, 



9 



we get the formula 



P A x(Ax;At) = J R JV (^) 



Ax | \ ~~ 3+7; — — 



x exp 



Ri 



AxW 2 -^f /|AxK^ 

+R2a v7m) 



'At 



0{a A \Ax\^=^ 



x 1 + R 



30 xm) 



2(2+a-2g) ' 
X 1 + AX 3 + a -2f, 



+ 



(33) 



O I a* \Ax\- *+°-'f> x 



where the constants Rjy, Ri, R2, -R3 contain the combination 4a / g 2 of the remaining 
parameters and depend also on a and /3. The most important aspect of this result is the 
stretched-exponential decay of the tail of the distribution. Therefore, it is fatter than 
the Gaussian, but still not fat enough to explain the power-law tails. 

There are several special cases to be investigated separately. First, it is clear that 
the algebraic substitution which leads to the function ((32)) fails ii (3 — 1/4. In this case, 
exponential substitution must be used instead, in the form 



ip{u;r]) 



.(a+3/2)« 



(a+l/2)« 



(34) 



which leads to a result very similar to (|33p . In fact, it is equivalent to simply setting 
P = 1/4 in @3). 

More fundamental changes occur in the cases we called previously "Heston-type" and 
"GARCH-type" . The saddle-point method is unnecessary here and the integral (|2T)]) can 
be performed explicitly. Thus, for 1 + a — 2/3 = we obtain 



P AX (Ax; At) 



2J\f 



(1-4/3 + Aaa/g 2 ) V2n At \ 
( \Ax\ [4a 



(Ax) 2 g 2 
At 4a 



(35) 



x Ki 



with K v (z) the modified Bessel function. From the asymptotic expansion of the Bessel 
function we deduce exponential decay of the tail of the return distribution. This result 
was already found for the Heston model itself 27, 28, l29j . Now we can see that this 
behaviour holds unchanged for the whole "Heston-type" class of models. 

For 2 + a — 2/3 = the calculation is even simpler. The result is the algebraic decay 
of the return distribution 



Pax (Ax; At) 



27V 



r( 



2/3 + 2a/ g 2 ) 



V2tt At 



(Ax) 2 
2At 



2a a 
9 2 
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-i+2/3+2a/ ff 2 



(36) 



valid for entire "GARCH-type" class. We have already seen that the moment method 
reveals the power-law tail at large times for the "GARCH" model, a = 0, /3 = 1. Now 
we can see that, at least for short time, the power-law tails apply to a wide ensemble of 
models. 



5.2.2. 7 = 1 

The same calculation can be performed for the case 7 = 1. For some values of a and 
(3 it is sensible to accept also negative values of the stochastic volatility St- This is the 
case, e.g. of the Stein-Stein model, a = (3 = 0. However, for generic values of a and 
(3 we still need to restrict the definition range to St > 0. Therefore, we postulate the 
reflecting boundary at St = also here. This makes the case 7 = 1 entirely analogous to 
the previously studied 7 = 1/2. 

The special cases are /3 = 0, which requires exponential substitution, and again the 
1 + a — 2/3 = and 2 + a — 2(3 — cases as before. The generic case is characterised by 
substitutions 



(Ax) 



77 = a 



At 

*g*_ (Ax) 2 
4a 



sr. 

4a 



(2 + a - 2(3) 
w (2 + 



3+ as— 2/3 
2/?)4+o-2£ 



G = 



1 



At 

g 2 (Ax) 2 
4a At 



2(3 



(2 + a - 2(3) 



20 



The integral to perform is 



where 



Pax(Az;A£) = 



AfG 
V2nAt 



cxp ( - -i){u; rjfj du 



up +u 2 " 
Then, the return distribution for Air — > 00 is 



rj u 



2P 



(37) 



(38) 
(39) 



P AX (Ax;At) =R N (\M 
\ v At 




(40) 



with, again, some constants Rn, Ri, R2, -R3, in full analogy with (|3"3"|) . 
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Explicit calculation shows that the case (3 — 0, treated by exponential substitution 
and saddle point method, gives again the same result as direct substitution /3 = in 
(US). 

The special cases are treated similarly. For 1 + a — 2/3 = we have the following 
substitutions 



G 



(Ax) 



At 



(-) 

\4aJ 



9 2 (Ax) 2 



4a At 



2£ , 2a 

" 3 3g 



-2(3 + 2aa/g 2 
and the function within the exponential is now 

IpOu) = u-Z/i-^+^/a 2 ) +u l/(-2l3+2aa/g 2 ) 



(41) 



(42) 



The resulting return distribution can be expressed by generalised hypergeometric func- 
tions, but it is more transparent to perform the saddle-point calculations, which give the 
tails 



P AX (Ax; At) =R 




exp 



(43) 



The remaining special case is 2 + a — 2(3 
form 

P AX ( Ax; At) 



0. The integral to perform has the following 
(Ax) 2 



l-2/3-2a/g^ 



exp i 



2aa 
g 2 s 



, . (44) 

x 2-AI J u ' ' V 2Ats 2 ' ' 

Here, the saddle-point method is of no use, because the corresponding function ip(u), 
if used, would have no local extreme. Instead, we can use a simple trick to estimate 
the behaviour of the tail, at |Ax| — > oo. The exponential function within the integrand 
is close to 1 for s > s* ~ \Ax\/y/At and very small for s < s* . More precisely, this 
behaviour holds only for \Ax\/y/At > 2aa/g 2 , but for the discussion of the tail it is quite 
sufficient. Thus, the integral (PHf can be estimated by integrating only the algebraic 
factor from s* to infinity. Finally, we obtain the power-law tail of the return distribution 
in the form 

(\Ax\\- 2 0~ 2a /9' 2 



Pax (Ax; At) 



(45) 



As a summary, we can see that the short-time return distribution has generically a 
tail in the form of a stretched exponential. This becomes simple exponential in the case of 
7 = 1/2 and Heston-type combination of parameters, l+a—2{3 = 0. The exceptional case 
is the family of "GARCH"-type models, characterised by the combination 2 + a — 2/3 — 0. 
In that case, we found that both for 7 = 1/2 and for 7=1 the tail of the return 
distribution has a power-law form. 

As a separate calculation, we show how the same procedure works for the expOU 
model. We have to perform the integral 



Pax (Ax; At) 



7T 2 At ff 2 
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(46) 



where 
and 



ip(u) = 2£u+ (In u)' 
9 2 (Ax) 2 



2a At 



(47) 



(48) 



The saddle point u* can be expressed using the Lambert function Wj J (y) defined by the 
equation Wl(j/) e Wh ^ — y. Finally, we find that 



^(«*) = (1 + WL(0V 
and within the saddle-point approximation 



g 2 (Ax) 



p«(^ao«-p{-^[(i + W l(S^L))-i]} 



(49) 



(50) 



The tails of the return distributions, Ax — > oo, correspond again to the limit £ — > oo. 
However, the situation is less transparent here, because the Lambert function does not 
have a simple asymptotic behaviour. Naively, one could use the asymptotics 

Wi,(y) — my, y->oo (51) 

but it holds well only for extremely large y. There is a much better approximate formula 

1 



W L (y) ~ lny 



1 + In y 



- 1 lnlny, 



y 



(52) 



which is reasonably accurate already for y > 3. The farthest tail of the return distribu- 
tion, using (|50l) and (|5"TT) . has the form 



Pax (Ax) ~ (Ax) 



-(2a/g 2 )lnAz 



(53) 



which is steeper than any power but fatter than any stretched exponential. Therefore, 
the expOU model lies somewhere between the "GARCH" model and the other models 
exhibiting stretched exponential tails. In this sense it is very promising, as the true 
power-law tails are hardly detectable in real data, but the behaviour very close to power 
law is well documented. The exponential Ornstein-Uhlenbeck process may serve as a 
convenient tool in this situation. 



6. Conclusions 

We showed that many properties of wide class of stochastic volatility models can be 
deduced quite generally without full explicit solution of the Fokker-Planck equation for 
the set of coupled stochastic differential equations. The class characterised by algebraic 
coefficients was solved systematically in the long-time and short-time regimes. Full so- 
lution for the stationary autocorrelation functions was found in two special cases, while 
generic features of the autocorrelation function are found generally. 

Specifically, the method of moments was shown to reveal the asymptotic distribution 
of return and volatility. The return distribution is asymptotically Gaussian, except for 
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the "GARCH" model, where power-law tail develops. This tail in returns originates from 
analogous power-law tail in the distribution of volatility. The empirical findings 471 148 1 
suggest that the central part of the volatility distribution is close to log-normal, while 
the tail is fatter, perhaps indeed a power-law. This makes the "GARCH" stochastic 
volatility model a good candidate for phenomenological modelling and risk assessment. 
Saying "phenomenological" we admit and stress that the model does not say anything 
of the real mechanism how the price fluctuations are produced, neither it explains why 
just the chosen combination of the parameters a, (3 and 7 is superior to others. 

While the price process is not Markovian any more, the compound price-volatility 
process does have the Markov property by definition. It results in exponential decay of 
volatility autocorrelation. Multiple time scales found in the expOU model may be used as 
explanation for the slow decay of empirical autocorrelations [H|,[l(|. However, to balance 
the opinions, we should also note that the approach of [l(| was severely criticised in [52j ]. 

The fat tails in return distribution, which is the most notorious stylised fact of em- 
pirical econophysics, are in fact a short-time property. At long time distances, the return 
distribution approaches to a Gaussian, as if the price process was a Brownian motion. 
This lead us to the study of short-time properties of the stochastic-volatility models, 
which at long times exhibit undifferentiated (with exception of the "GARCH" class) 
Gaussian behaviour. We examined the return distributions for general values of the pa- 
rameters a, (3 and 7 and we found that stretched exponential is the generic form of the 
distribution. As a special case, we find exponential distribution in the Heston-type class, 
1 + a - 2(3 = 0. One exception is the "GARCH" -type class, (2 + a - 2(3 = 0) where 
the distribution has power-law tail. Of course, this is not surprising, as we already know 
that the power-law tail persists for all times. The intermediate class is represented by 
the exponential Ornstein-Uhlenbeck process, where the tail is thinner than a power-law, 
but fatter than stretched exponential. 

To sum up, we examined systematically whole class of stochastic volatility models 
and found explicitly their asymptotic and short time properties, as well as the volatility 
autocorrelation. Power-law tails are reproduces in the family of "GARCH" -type models, 
while generic form of the short-time return distribution is stretched exponential. 
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